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In this paper we analyze a coupling between the very large jumps in physical and operational times 
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superdiffusive features of anomalous diffusion underlying the two-power-law relaxation patterns. 
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I. INTRODUCTION 

The Continuous Time Random Walk (CTRW) formal- 
ism is a very powerful stochastic approach to model phys- 
ical processes demonstrating anomalous diffusion and 
slow, power-law, relaxation. It describes random walks in 
space and time by means of iid (independent and identi- 
cally distributed) couples of space and time random steps 
(Ri,Ti). The simplest, decoupled CTRW considers in- 
dependent time and space steps. This model involves 
stable distributions, and it shows various anomalous be- 
haviors like subdiffusion (diffusion slower than normal 
one), Mittag-Leffler relaxation and fractional diffusive 
equations [l|-|3| . A more complex CTRW model accounts 
for coupling between time and space steps. The coupled 
CTRWs were considered in the context of anomalous dif- 
fusion and non-exponential relaxation • In this case 
the anomalous diffusion evolution is much richer. Sub- 
and superdiffusion (faster than normal) can be modeled. 
However, the analysis is rather exotic for the research, 
and it is in progress. Recently, the anomalous subdiffu- 
sive behavior attracts a great attention in modeling of 
subdiffusion in space-time-dependent force fields beyond 
the fractional Fokker-Planck equation 0, Q. This ap- 
proach uses the Langevin-type dynamics with subordina- 
tion techniques, where the force depends on a compound 
subordinator. It is coupled because of a Levy-stable pro- 
cess directed by its inverse. The fractional two-power- 
law relaxation can be also described in the framework of 
coupled CTRWs based on subordination of a stochastic 
process with the heavy-tailed distribution of the waiting 
times by its inverse [I [Tol. Although the papers have a 
different physical background, they intersect into the ap- 
plication of the coupling between the Levy-stable process 
and its inverse. Undoubtedly, this new random process 



(subordinator) is of an essential interest for understand- 
ing of the anomalous relaxation phenomena and was in- 
vestigated unsufficiently yet. In this paper we are going 
to make up the deficiency. 



II. COUPLING BETWEEN THE VERY LARGE 
JUMPS IN PHYSICAL AND OPERATIONAL 
TIMES 

The probability density of the position vector Yt = Bs t 
(where B T is the standard Brownian motion) can be 
found from a weighted integration of the joint probabil- 
ity density of the couple (R r ,T r ) over the internal time 
parameter r by subordination. The stochastic time evo- 
lution T T and its (left) inverse process St permits one to 
underestimate or overestimate the physical time t. 

The sum of iid heavy-tailed random variables Tj 



Pr(2i>*)~ ± 



as t 



(1) 



< a < 1, to > converges to a stable random vari- 
able in distribution as the number of summands tends to 
infinity. Let U n = X)"=o w ith — 0. The counting 
process N t = max{n G N | U n < t} is inverse to U n which 
can be defined equivalently as the process satisfying 



U Ni < t < U Nt+1 for t > , 



(2) 
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what follows directly from its definition. In fact, the two 
processes f/jy t and C/jv t +i correspond to underestimating 
and overestimating the real time t from the random time 
steps Ti of the CTRWs. 

In terminology of the Feller's book [ll| the variable 
Z t = UN t +i — t is the residual waiting time (life-time) 
at the epoch t, and Y t = t — Ujf t is the spent wait- 
ing time (age of the object that is alive at time t). 
The importance of these variables can be explained by 
one remarkable property. For t — > 00 the variables Y t 



2 



and Z t have a common proper limit distribution only 
if their probability distributions F(y) and F(z) have 
finite expectations. However, if the distribution Fix) 
satisfies 1 — F(x) = x~ a L(x), where < a < 1 and 
L(xt)/L(t) — > 1 as x — > oo, then according to the 
probability density function (pdf ) of the normalized vari- 
able Yt/t is given by the generalized arc sine law 



sin(Tra) 1 

Pa\x) = X (1 — X) 



while Z t /t obeys 



, , sin(7ra) „ . . -, 



(3) 



(4) 



Since Ejv t = t — Y t and Sjv t +i = Z t + 1, the distributions 
of Ejv t /i and Sjv t +i/i can be obtained from Eqs. © 
and ^ by a simple change of variables 1 — x = y and 
1 + x = z, respectively. 

We now return to the processes J7jv t and UN t +i in- 
troduced above. Recall that Ti are iid positive random 
variables with a long-tailed distribution ([l) . In this case 

UN t /t tends in distribution (— >) in the long-time limit to 
random variable Y with density 

p Y ( x ) = S j?(™)- x <*-\l- x )-°<, <X <1 (5) 



and UN t +i/t —> Z with the pdf equal to 



x > 1. 



(6) 



| T r (Ar) > t} or equivalently S t (Ar) = At N t . There- 
fore, in the limit At — > the processes E/jv t +i an d i7jv t 
can be expressed through the stochastic process T T and 
its left limit subordinated by their inverse: 

U Nt -4 T St and U Nt+ i A T St . 

The passage from the discrete process Tj to the continu- 
ous one T T allows one to reformulate the unequality © 

as 



T St <t< T St for t > . , 



(7) 



underestimating or overestimating the real time t. From 
Theorem 1.13 in [l3[ the joint probability density p(y, z) 
of T St and T St with < T<7 < t < T St takes the form 

asin(7ra) „ -,. , i „ 

for < y < t < z. After integrating © with respect to 
z in the limits [t , oo[ (or with respect to y in the limits 
[0 we obtain the densities of Tg and Ts t , respectively 



V (t,y) 

p+(t,z) 



sin na 



sm not 



v a -\t-y)- a , o< y <t, (9) 
z- 1 t a { z -ty a , z>t, (io) 



valid for any time t > (see Fig. [TJ . The moments of Tg 
and Tg t can be calculated directly from the moments of 
Y and Z by using relations 



The functions p Y (x) and p z (x) correspond to special 
cases of the well-known beta density. It should be no- 
ticed that the density p Y (x) concentrates near and 1, 
whereas p z {x) does near 1. Near 1 both tend to infinity. 
This means that in the long-time limit the most probable 
values for Ujf t occur near and 1, while for Upf t +i they 
tend to be situated near 1. 

As a consequence, the random variable Y has finite 
moments of any order. They can be calculated directly 
from the density (J5J) and take the form 

(Y) = a, (Y 2 ) = — ^— ^— — ^ , • • • j 



a(l + a) . . . (a + n — 1) 



where n£N, while even the first moment of Z diverges. 
The divergence of t/jv t +i results from the long-tail prop- 
erty of the time steps Tj so that (Tj) = oo, yielding 
too long overshot above t. 

The nonequality @ can also be represented in a 
schematic picture of time steps as T~(At) = U[ t /at] an d 
T t (At) = U[ t /at]+i, where T r "(AT) = lim^o T T _ e (AT) 
is the left-limit process, and [x] indicates the integer part 
of the real number x so that [x] < x < [x] + 1. The in- 
verse process of T~( At) andT r (AT) is St (At) = inf{T > 



Tg t = tY and T St = tZ . 



where = means the equality in distribution. Thus, the 
process Tg has finite moments of any order, while Tg t 
gives us even no finite the first moment. The overshot of 
Ts t > t is too long also in the limit formulation. Notice 
that p + (t,y) = y~ 2 p~ (t^ 1 ,y~ 1 ). At this point we should 
mention that compound subordinators, and in particular 
the subordination by an inverse Levy-stable process via 
a Levy-stable process, were considered already in [l4l |. 
However, the construction of compound subordinators 
has been based on the statistically independent stochas- 
tic processes. This leads to quite different results in com- 
parison with ours. In our construction of the compound 
subordinators Tg and Ts t the processes U t and S t are 
clearly coupled. 



III. ANOMALOUS DIFFUSION WITH UNDER- 
AND OVERSHOOTING SUBORDINATION 

According to [10(, the widely observed fractional two- 
power relaxation dependencies 



for u»(j„ 



(11) 
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FIG. 1: (Color online) The probability density p~(y) with 
support on < y < t and the density p + (y) with support on 
y > t for different values of the index a. 

and 

Ax(w) ~ (iuj/ujp)™ 1 for ui <C w p (12) 

of the complex susceptibility x( w ) — — *x"( w )j 

where A%(w) = x(0) — x(^), the exponent n and m fall in 
the range (0, 1), and uj p denotes the loss peak frequency, 
are closely connected with the under- and overshooting 
subordination 

Z£ 7 (t) < S a (t) < Z£ 7 (t) for i>0, 

where Z£ 7 (i) = X? [£«,(*)], Z°„(t) = X? [$*(*)]• Here 
the processes X^{t) and -X^i) are nothing else as T<7 
and Tg t with the index 7. They are subordinated by 
an independent inverse a-stable process S a (t) forming 
the compound subordinators (t) and Z° (t), respec- 
tively. The approach enlarges the class of diffusive sce- 
narios in the framework of the CTRWs. This new type 
of coupled CTRWs follows from the clustering-jump ran- 
dom walks idea [HI]. As it has been rigorously proved 
[l6|, the clustering with finite- mean- value cluster sizes 
leads to the classical decoupled CTRW models, but as- 
suming a heavy-tailed cluster-size distribution with the 
tail exponent < 7 < 1, the coupling between jumps and 
intcrjump times tends to the compound operational times 
Z^ (i) and Z°(t) as under- and overshooting subordi- 
nators, respectively. 

The overshooting subordinator yields the anomalous 
diffusion scenario leading to the well-known Havriliak- 
Negami relaxation pattern [17], and the undershooting 
subordinator leads to a new relaxation law give n by 
the generalized Mittag-Leffler relaxation function [9|, [lfj . 
These results are in agreement with the idea of a super- 
position of the classical (exponential) Debye relaxations. 
Thus, the stochastic mechanism underlying the anoma- 
lous relaxation is quite clear, but the corresponding diffu- 
sion analysis requires some additional clarity. Let B(t) be 



the parent process that is subordinated either by Z^ (t) 
or Z®(t). Then the subordination relation, expressed 
by means of a mixture of pdf 's, takes the form 

/>oo />oo 

p r {x,t) = / / p B (x,y)p ± (y,T)p s (T,t)dydT , 
Jo Jo 

(13) 

where p r (x,t) is the probability density of the subor- 
dinated process B[Z%^(t)] (or B[Z^(t)\) with respect 
to the coordinate x and time t, p B (x,r) the probabil- 
ity density of the parent process, p ± (y,r) the probabil- 
ity density of T<7 and Ts t respectively, and p (r, t) the 
probability density of S(t). Recall that for the subdiffu- 
sion B[S(t)], by taking the Laplace transform from the 
corresponding subordination relation, we can derive the 
celebrated fractional Fokker-Planck equation [Tsj ] . It is 
therefore reasonable to ask is it possible to find a dif- 
fusion equation corresponding to relation (1131) . In the 
Laplace space 

/>oo 

/(«)= / e- ut f(t)dt 
Jo 

we obtain 

r°° dr 

p r (x,u)=u<*- 1 p B (x,u a /z)p+(z)- (14) 
Ji z 

with Pq(z) — sin(7T7) z _1 (z — 1) 1 /it for z > 1, as well 

p r (x,u) = u a - 1 f 1 p B (x,u a /z) Po -(z)- (15) 
Jo z 

with Pq(z) — sin(7T7) z 7_1 (l — z)~ 7 /n for < z < 1. 
The Laplace image of the pdf of the subordinated process 
B[S(t)] can be simply expressed in terms of an algebraic 
form with the Laplace image of the parent process pdf. 
This allows one to get the fractional Fokker-Planck equa- 
tion driving the spatio-temperal evolution of the propa- 
gator of the anomalous diffusion underlying the Mittag- 
Leffler relaxation 0, 0, Hq . However, expressions (fit)) 
and (|15l) are not similar to the latter. They have an inte- 
gral form. Nevertheless, derivation of the corresponding 
Fokker-Planck equation is also possible. 

If we take the Laplace transform with respect to t and 
the Fourier transform with respect to x for p r {x, t) in Eq. 
(fl3l) . the Fourier-Laplace (FL) image reads 

FL (p r )(k,s) 

= s"- 1 / / e~^ k)y p ± (y,T)e- Tsa dydr, (16) 
Jo Jo 

where tp(k) is the log-Fourier transform of the parent 
process pdf p B (x, y). Consider the case of p~{y, r). After 
changing variables y — zt we take the integral 

r°° 1 

/ e -r(s^(k) Z ) dT = . 

Jo s a +i;(k)z 
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Next, the change of variables t — z/(l — z) maps [0,1] 
onto [0 , oo). This helps to derive 



FL (p r )(M 



r(7)r(i-7)y (s a + ip(k))t + s a 



The last expression can be easily calculated from the in- 
tegral Q 



ti- 



■ dt = T(rf)r(i-i). 



io t + 1 

The FL image of p r (x,t) with the undershooting direct- 

'ly(t)=X!> 

a-y-1 



ing process Z^ At) — X^[S a (t)] is of the form 



FL (p r )(k,s) 



(s a + ip(k)) 



7 



(17) 



Finally, we invert the Fourier and Lapace transforms to 
get the pseudo-differential equation 



L FP (x) 



p r (x,t)=S(x) 



t -art 



r(l -ay) ' 



(18) 



where L FF (x) is the Fokker-Planck operator, S(x) the 
Dirac function, and d a /dt a denotes the Riemann- 
Louiville derivative. The corresponding Fokker-Planck 
equation can be obtained also in the case when the over- 
shooting directing process Z°(t) — X°[S a (t)] is taken 
into account. Unfortunately, the derivation is more com- 
plicated as we present below. 

In the case of p + {y, r), after the substitution y — zr, 
we map [1 , oo) onto [0 , 1] by the change of variables 
z = 1/x. Then we obtain the corresponding FL image 



FL (p r )(k, S ) = 



r(7)r(i - 7) Jo 



1 xf' 1 (1 - x) _ T dt 
s a + ip(k)/x 



The mapping t = x/(l — x) transforms the latter expres- 
sion to the form 



FL (p r )(k,s) 

„a-l 



V dt 



r( 7 )r(i - 7 ) J (i + w + m)t + 1>(k)] 

This integral can be calculated exactly: 



■dt = 



r( 7 )r(i- 7 ) 



1 - ib/aV 



l (t + l)(at + b) (a-b) 
As a result, the FL image of p r (x,t) with the directing 



process Z° (t) = X°[S a (t)], can be written as 



FL (p r )(k,s) 



1 



1 



i>(k) 



7 



(19) 



Now we invert the Fourier and Laplace transforms to get 
the pseudo-differential equation 



dt^ 



+ L FP (x) 



p r (x,t) = f a Jx,t), (20) 






7 






Lfp(x) 




L FP (x) 


') 



■ S(x) 



is a function depending on the probability density 
p B (x,y). The exact form of f an (x,t) is quite different 
from the right-side term of Eq . (TTgl) . In this connection it 
should be pointed out the work [20| , where the derivation 
of a fractional Fokker-Planck underlying the Havriliak- 
Negami type of relaxation is based on the entirely phe- 
nomenological approach of [2lJ. However, the stochastic 
background leading to the anomalous diffusion yielding 
the Havriliak-Negami pattern, has remained behind these 
works. It should be noticed that Eqs. ([T51) and (|2"0)) have 
been derived independently in papers [22| . |23| . 

To calculate the moments of the processes B[Z^ (t)] 
and B[Z° (t)], assume for simplicity, that the parent 
process B is a one-dimensional Brownian motion. Its 
moments are written as 



hn(t) 



1 



V^nDt 
(2n)! 



x 2n exp 



x 



dx 



(Dt) n , 

where D is the diffusion coefficient. If the subordinator 
Zyt) governs the Brownian motion, then the moment 
integral reads 



< x 2n > 



x 2n p r (x,t)dx 

= B n f z n PQ -(z)dz f 
Jo Jo 

= (2n)I Jjn (7,n) t no 



_n 



p (r,y)dr 



n\ T(l + not) ' 



(21) 



where (7, n) = 7(7 + 1)(7 + 2) . . . (7 + n — 1) is the Ap- 
pell's symbol with (7,0) = 1. When another subordina- 
tor Z® At) is used, even the first moment of the subordi- 
nated process B[Z° (t)] diverges because the probability 
density p$ (z) gives no finite moments. Thus, the process 
B[Z^ 7 (t)} is subdiffusion, and B[Z°^{t)} is superdiffu- 
sion. In Fig. [5J as an example, the propagator p r (x,t) 
for the under- and overshooting anomalous diffusion with 
a — 2/3 and 7 = 2/3 is drawn. 

It should be noticed that the ordinary subdiffusion 
B[S a (t)] takes an intermediate place between the under- 
and overshooting anomalous diffusion B[Z^ 7 (t)] and 
B[Z°„(t)]. The feature is illustrated in Fig. ED This al- 
lows one to compare an asymptotic behavior of the tem- 
poral evolution of diffusion fronts. From that one can see 
that the diffusion front of B[Z% (t)] is more stretched 
than the front of B[S a (t)), whereas the diffusion front 
of i3[Z^ 7 (t)] is more contracted in comparison with the 
front oi'B[S a {t)]. 

One of interesting questions is what interpretation can 
be assigned to the subordinators Z% (t) = X^[S a (t)] 
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FIG. 2: (Color online) Propagator p r (x,t) of the under- and 
overshooting anomalous diffusion with a constant potential, 
a = 2/3 and 7 = 2/3, drawn for consecutive dimensionless 
instances of time t — 1,3, 10. The cusp shape of the pdfs 
appears. 



ered separately. The inverse Levy-stable process S a (t) 
accounts for the amount of time, when a walker does not 
participate in motion. The pdf of the subordinated pro- 
cess B[XjJ (t)] is a special case of the Dirichlet average, 
namely 

F(j,x,t) = [ P b {x,tz)z^ 1 (1- z)" 1 dz. 

n Jo 

Recall that many of important special and elementary 
functions can be represented as Dirichlet averages of con- 
tinuous functions (see more details in 24] ) . The Dirichlet 
average includes the well-known means (arithmetic, geo- 
metric and others) as special cases. The process X^(t) 
evolves to infinity like time t. Its contribution in the sub- 
ordinated process B[X^ (t)] is taken into account by the 
Dirichlet average of the probability density of the parent 
process B. The similar reasoning can be developed for 
the process X^(i). 



IV. CONCLUSIONS 




-6 -4 -2 2 4 6 



X 



FIG. 3: (Color online) Left panel: the propagator p r (x,t) 
of under- (a) and overshooting (b) anomalous diffusion with 
a — 2/3 and 7 = 2/3 for t = 1. The line (c) displays the 
propagator of ordinary subdiffusion with a — 2/3 and 7 — 
1 for t — 1. Right panel: diagram shows the interrelation 
between B[Z„ n (t)], B[S a (t)] and B[Zg n (t)]. Here m and 
l — Ti denote the power-law exponents defined in formulas 
(fTT]) and gll). 

and Zg 7 (t) = X°[S a (t)]. As the processes X^(r) and 
X^(t) are independent on S a (t), they can be consid- 



The paper introduces an approach to study of the cou- 
pling between the very large jumps in physical and opera- 
tional times. It is based on the compound subordination 
of a Levy-stable process T(t) by its inverse S(t). The 
inverse Levy-stable process is actually the left-inverse 
process of the Levy-stable process. In fact, we have 
S[T(t)} = r, while T[S(t)] > t holds. In the framework 
of CTRWs and the Langevin-type stochastic differential 
equations the compound subordinator provides a direct 
coupling of physical and operational times. The subor- 
dination scenario leads to two types of operational time: 
the spent life-time and the residual age. In the first ran- 
dom process all the moments are finite, whereas the sec- 
ond process has no finite moments. We have shown that 
the approach is useful for analysis of anomalous diffusion 
underlying all empirical fractional two-power-law relax- 
ation responses. Due to the two types of the operational 
time the diffusion can display as well the subdiffusive and 
superdiffusive character. 
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